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Abstract 

The radiative B — > t^t^'j decay is investigated in the framework of the Standard 
Model. When only short (short and long together) distance contributions are taken 
into account, the Branching Ratio is found as 9.54 x 10^^ (1.52 x 10^^), for the value 
of the cut 6 = 0.01 imposed on the photon energy. 



1 Introduction 



Experimental discovery of the inclusive and exclusive B — > Xs'j and B K*^ |l|] decays 
stimulated the study of the radiative rare B meson decays with a new momentum. From 
theoretical point of view this is due to the fact that they are very sensitive to the flavor 
structure of the electroweak interactions, as well as QCD radiative corrections and the new 
physics beyond the Standard Model (SM) [0]. From experimental point of view studying 
radiative 5-meson decays allows more precise determination of the parameters of the SM, 
such as the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix, the leptonic 
decay constants etc., which are yet poorly known (see for example 0). 

Currently there is impressive effort in this direction, and many new facilities are under 
construction for studying the rare 5-meson decays, namely, symmetric and asymmetric 
S-meson factories at Cornell, KEK and SLAC. Progress is also being made in hadronic 
environment at HERA-B and there are some plans for TeV-B and LHC-B. These ma- 
chines will serve to measure the processes, for which SM predicts very small Branching 
Ratios. Among the rare decays, the flavor changing decays of the 5-meson which proceed 
via electroweak penguins, are of special interest due to their relative cleanliness and their 
sensitivity to the new physics. The rare B t^t~'^ decay belongs to this category. 

From helicity arguments it is clear that the matrix element ol B i^i^ {i = e, /i) 
decay will be proportional to the lepton mass and therefore the corresponding Branching 
Ratios will be strongly suppressed. Note that in SM, B{B — > e~^e~) ~ 4.2 x 10~^^ and 
B{B — > ~ 1.8 X 10"^. It is well known however that the the B{B —>■ r+r^) ^ 8 x 10"'^ 

in SM [Q, and thus this decay can be measured in future i?-factories with high enough 
efficiency. 

When photon is emitted in addition to the lepton pair, no helicity suppression exists, 
and a "large" Branching Ratio is expected. Indeed in it was shown that the B{B — > 
e"'"e~7) ~ 2.35 x 10~^. For B {i = e, /i), the contributions of the diagrams, 

where photon is radiated from any charged internal line, can safely be neglected, as they 
are strongly suppressed by a factor m^/m^^/ in the Wilson coefficients (see 0). Moreover, 
it follows from helicity arguments that, the contributions of the diagrams where a photon 
is emitted from the final charged lepton lines, must be proportional to the lepton mass 
rrii {i = e, /x), and hence they are neghgible as well. Therefore in B {£ = e, fi), 

the main contribution should come from the diagrams, where photon is emitted from the 
initial quarks. 

In 5 — > r"'"r~7 decay, the situation becomes very different. In this case, we cannot 
neglect the contribution of the diagrams, where photon is radiated from the final r-leptons, 
since the mass of the r-lepton is not so much smaller than that of the 5-meson. So, in 
B T^T~'-) decay comparable contributions come from diagrams where photon is radiated 
both from initial and final fermions. These contributions can give essential information 
about the relative roles of the strong and electroweak interactions. 

In this work we investigate the B — > r"'"r~7 decay, and the paper is organized as follows. 
In Section 2 we give the necessary theoretical framework for the B — > r+r~7 decay. Section 
3 is devoted to the numerical analysis and the discussion of the results. In Appendix the 
detailed description of the cancellation of the infrared (IR) singularities, is given. 
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2 Theoretical framework for the B t^t 7 decay. 

The matrix element for the b sr^r^'-f decay can be obtained from that of the b —>■ st^t^ . 
It is well known that the short distance contributions to 6 ^ sr+r^ decay comes from the 
box, Z- and photon-mediated penguins. Thus, in the SM, QCD-corrected amplitude for 
b ST^T~ can be written as 0, 0]- 



M 



2V2 



f:^ytbVtli^Cg^^s-fi,{l - 75)6 r7^r + CioS7/,(l - 75)&^7/.75r 



- 2C7^sia^^pi,{l + 75)&r7^r | . (1) 

In Eq. (1) p is the momentum transfer, and the mass of the strange quark is neglected, 
Vi/s are the corresponding elements of the CKM matrix. The analytical expression of all 
Wilson coefficients Cg^^ , Cio and C7 can be found in @, 0. In order to obtain the matrix 
element for b sr^r^'y, it is necessary to attach photon to any charged internal, as well 
as external line. Contributions of the diagrams with photon attached to the any charged 
internal line, are strongly suppressed and therefore we shall neglect these in the following 
discussions. Thus, as explained previously the main contributions to the b sT^T~'~f decay 
comes from diagrams, when photon is radiated from initial and final fermions. 

When a photon is attached to the initial quark lines, the corresponding matrix element 
for the B — > t^t'-^ decay can be written as 



Ml = {i\M\B) 



+ C'iot'7m75^(7|s7m(1 - l5)b\B) 

- 2C7^(7|sza^.p.(l + 75)&|i?)T7Mr| • (2) 

These matrix elements can be written in terms of the two independent, gauge invariant, 
parity conserving and parity violating form factors |^ H: 



{j\sj,{l - j,)b\B) 




+ z 



- fiP ) 



+ i 



<{pq)-{e*p)q, h{p')\ . (3) 
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Here and are the four vector polarization and four momentum of the photon, respec- 
tively, and p is the momentum transfer. Substituting Eq. (3) in (4), for the matrix element 
M.I (structure dependent part) we get 



Ml = 



2^2 TT 



+ i 



(4) 



where 



A 
B 
C 
D 



Cl"g{p')-2C,'^g,{p^) 

ci^'m-2c,^Mp') 



ClO ^ f 2 



ml 



C 



10 



nil 



9{P ) , 
/(/) ■ 



(5) 



When a photon is radiated from the final r-leptons, the corresponding matrix element 
is (Bremstrahlung part) 



M2 = ^ r- VtbVtleifBCio2mr 

ZV 2 TT 



\ 2piq 2p2q J ^ 



(6) 



where Pb is the momentum of the S-meson. In obtaining this expression we have used 

(0|S7^75&|5) = -tfBPB,., 
(0|sa^,(l+75)6|5) = 0, (7) 

and the conservation of the vector current. 

Finally, the total matrix element for the B t~^t~j decay is obtained as a sum of the 
Ml and M2- 



M^ M1 + M2 . 



(8) 



The square of the matrix element, summed over the spins of the r-leptons and the polar- 
ization of the photon, can be written as 



\M\^ =\Mi\^ + \M2\^ + 2Re{MiM;) , 



(9) 



where 
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aGf 



2^2 TT 

+ 4 
+ 4 



Cf+\Df 



A\^ + \B\^ 



Ana js Re {B*C + A*D) {p^q - p2q) {piq + P2q) 

(f - 2ml) iiPKlf + (P2?)') - 4m^ (pig) {p^q) 
p" + 2m^) {{p^qf + (p2?)') 



+ Ami (PiQ) (p^q) 



(10) 



2Re(MiM;) 



aGr 



2V2 



71 



Ana I 16 Ciofsml 



Re{A) 



+ Re(D) 



(Piq) {P2q) 



{piq + P2qf 
(PiQ) {P2q) 



(11) 



aGf 



+ 



+ 



2V2n 
16 r 

Piq 

16 



VtbVZ 



Anal -32 + 



(Pi?)' 



3m^ -p^ - P2q 



+ 8-^{p' + 2p,q) 



P2q 



3ml -P^ - Piq 



+ 8 



p^ 



(Piq) {P2q) 



2ml - P 



(12) 



Here pi, p2 are momenta of the final r-leptons, and q is the photon momentum. The 
quantity depends only on the scalar products of the momenta of the external particles. 
In the rest frame of the 5-meson, all these scalar products are fixed, if the photon energy 
and the lepton energy Ei are specified. The Dalitz boundary is given as 



0<E^< 



m 



R - Ami 



2mB 



(13) 



7Tij5 — E^ E^ 



1 - 



Ami 



2 '\ m% — 2mBE^ 



<Ei< 



ms — E^ 



E^ 



1 - 



Ami 



m^ — 2mBE^ 



(14) 



The \M.i\ term is completely infrared-frcc; the interference term has an integrable infrared 
singularity. Thus infrared divergence appears only in |7Vf2|^- The infrared singularity 
originates in the Bremstrahlung processes, when photon is soft. It is clear that in this 
hmit, the B — > t+t~7 decay cannot be distinguished from B — > r+r". Therefore both 
processes must be considered together in order to obtain finite result for the decay rate. 
In the Appendix we show that IR singular terms in |A^2|^ exactly cancel the 0{a.) virtual 
correction in i? ^ r+r^ amplitude. 

In this work, our point of view is slightly different from the standard description. 
Namely, we consider the Bremstrahlung process as a different process but not as the a 



A 



correction to the B — >■ r+r" decay. In other words, we consider the photon in i? ^ t^t^'j 
as a hard photon. Therefore, in order to obtain the decay width of the B — > r^r^ + 
(hard photon), we must impose a cut on the photon energy, which will correspond to the 
experimental cut imposed on the minimum energy for detectable photon. We require the 
energy of the photon to be larger than 50 MeV, i.e., > arriB, where a > 0.01. 

After integrating over the phase space, and taking into account the cut for the photon 
energy, for the decay rate we get. 



2V2 



TC 



a 



(2 7r) 



moTT 



3 '"B 



X 



1 

12 Js 

2 



l-4r 



X dx \ 1 



4r 



\D\ 



V 1 — X 

[1- X - 4r) 



m 



B 



Af + \Bn (l-x + 2r) 



l-4r 



1 + wi- 



- 2CiofBr / x^ dx Re (A) In- 



4r 
1 — X 



4r 
1 — X 



4|/ijCio|'r^ 



l-4r 



dx 



4r 2 
2H x] In 

X X 



1 + Wl 



4r 
1 — X 



4r 
1 — X 



+ -{l-x) 

X 



4r 
1 — X 



where x 



2E^ 



— - is the dimensionless photon energy, r 
niB 

4:ml 

S<x<l ^ . 



mz 



m 



and 5 = 2a, satisfying 



B 



m 



B 



(15) 



From Eq. (15) it follows that for calculating the decay width we need explicit forms of 
the form factors g, /, gi and /i. These form factors are calculated in the framework of 
light-cone QCD sum rules in |Q (see also [§), and their dependences, to a very good 
accuracy, can be represented in the following dipole forms. 



9{P 



1 GeV 



:i- 



„2 ' 

5.62^ 



0.8 GeV 



(1- 



_P \2 

6.52^ 



, 2x 3.74 GeV^ , , 0.68 GeV^ , , 

(1_^)2 (1--)' 
^ 40.5^ ^ 30^ 

which we will use in the numerical analysis. 
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3 Numerical analysis and discussion 

For the input parameters, which enter into the expression for the decay width we use the fol- 
lowing values: nib = 4.8 GeV, nic = 1.35 GeV, = 1.78 GeV, tub = 5.28 GeV, \VtbVts* \ = 
0.045. We use the pole form of the form factors given in Eq. (16). For 5-meson life time we 
take r (Bs) = 1.64 x 10~^^ s Q. The value of the Wilson coefficients Gj {nib) and Cio (nib), 
to the leading logarithmic approximation, are (see for example [§], 0)- 

Gj (nib) = - 0.315, Gw (nib) = - 4.6242 . (17) 

The expression Gg^-^ (nib) for the b ^ s transition, in the next-to-leading-order approx- 
imation, is given as 

G^^^inib) = Gg{mb) + 0.l24w{s)+g{m„s){3Gi + G2 + 3Gs + G4 + 3G5 + Ge) 
-\g{ni,, s)(C3 + 3C4) - ]^g{rhb, s)(4C3 + 4^4 + 3^5 + G^) 

+ ^(3C3 + C4 + 3C5 + C6) , (18) 

with 

Ci = - 0.249, C2 = 1.108, C3 = 1.112 X 10"^C4 = - 2.569 X 10~^ 

= 7.4 X 10~^ Gq = - 3.144 X 10"^ Cg = 4.227 , (19) 

where rhq = mg/mb, s = /ml. The value of G^^^ for the h ^ d transition can be obtained 
by adding to Eq. (18) the term 

A„ [g {rhc, s) - g {rhd, s)] (3Ci + G2) , 

where 



A. 
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and replacing VtbV^^ in Eq. (15) by VtbV^^. In Eq. (18), w{s) represents the one gluon 
correction to the matrix element of operator O9, and its explicit form can be found in 
1^, while the function g{niq, s) arises from the one-loop contributions of the four quark 
operators Oi-Og, i.e.. 



g{mg, s') = --ln(mg) + — + -yg - -{2 + yg)^\l - yg\ 



X I 0(1 - yg) (\n ]^'^Ljl - m ) + e{yg - l)arctg- 



where yg = 4m^/s', and s' = /ml. 

For a more complete analysis of the B r^r^'y decay, one has to take into account the 
long distance contributions. For this aim it is necessary to make the following replacement 



g{ni,, s') ^ g{m,, s') - — ^ \, . . L'"' , (20) 
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where my = mv/rrif,, Ttot = ^/fTib- 

Our results for the Branching Ratio B{B — > t~^t~j) for two different values of the cut 
{S — 0.01 and S — 0.02) are presented in Table 1. 





Short Distance 
Contributions to the 
Branching Ratio 


Short and Long Distance 
Contributions to the 
Branching Ratio 




S = 0.01 


5 = 0.02 


5 = 0.01 


5 = 0.02 


Structure 


4.19 X 10-9 


4.19 X 10-9 


9.95 X 10-9 


7.68 X 10-9 


Bremstrahlung 
part 


4.11 X 10-9 


3.16 X 10-9 


4.11 X 10-9 


3.16 X 10-9 


Interference 
part 


1.24 X 10-9 


1.23 X 10-9 


1.16 X 10-9 


1.16 X 10-9 


Total 


9.54 X 10-9 


8.59 X 10-9 


1.52 X 10-s 


1.20 X 10-s 



Table 1: 



Note that, when only short distance effects are taken into consideration, the structure- 
dependent and the Bremstrahlung parts give, more or less, comparable contributions. How- 
ever, if long distance contributions are also considered in addition to the short ones, the 
structure-dependent part contributes more as compared to that of the Bremstrahlung part. 
This is due to the fact that, the structure dependent part contains {J/ip, resonance con- 
tributions in the latter case (see Eq. (20)). 

Prom these results it follows that there is a good chance for detecting r lepton decay in 
the future S-meson factories, provided that the efficiency is ~ 1/3. 
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Figure Captions 

1. Dependence of the Differential Branching Ratio for the B — > t'*"t~7 decay on the dimen- 

2E 

sionless variable x — — -, for the value of the cut 5 — 0.01 imposed on the photon energy. 

rriB 

In this figure, the curve with the sharp peak represents the long distance contributions. 

2. Same as Fig. 1, but for the cut value 6 — 0.02. 
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Fig. 1 
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Fig. 2 
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Appendix : Cancellation of the infrared divergences 



In this section we will show the explicit cancellation between the infrared singularities 
arising from the Bremstrahlung of a soft photon in the B —>■ t^t^'j rate and in the 0{a) 
virtual corrections to the B — >■ r+r" amplitude. 

As we have explained in Section 2, we are not going to include 0{a) virtual corrections to 
the B T^T^ amplitude in the calculation of the rate oi B ^ t^t^^. The Bremstrahlung 
process here is not considered as an 0{a) correction to the B t^t~' amplitude, but as a 
different process, namely the decay of a B-meson into r-lepton pair plus a hard photon. In 
order to calculate the physical rate of interest, we just have to impose a cut on the energy 
of photon, which will correspond to the experimental cut imposed on the minimum energy 
for detectable photon. 

In this sense our approach is very similar to the one given in [0, in studying Bg 
Xs'-ff decay. Therefore, in the present Appendix we will consider only those aspects of the 
discussion which are necessary to show the cancellation of the infrared singularities. For 



removing the infrared singularities we will use the method presented in , where the main 
idea is to use the dimensional regularization method, i.e., to replace 



2E^{2tiY 2|g|(27r) 



(n-l) 



(A.l) 



1/2 



where q = {q^, qt) is an n dimensional light-like vector: 

(f = kl = H ^ iQn- 

All calculations are performed in the rest frame of the 5-meson. We are interested in the 
situation when charged lepton is observed in the energy interval Em — Ai? < E < Em, 
where Em is that of the charged lepton energy in the two-body decay B t^t^ and 
/S.E « Em- We will retain terms of logarithmic and zeroth order in AE. It is clear that in 
this limit, only soft photons give contributions. 

Let us first consider the Bremstrahlung part. From |A^2|^ it follows that only terms 
proportional to 1/ (pig)^ , 1/ (P2Q')^ and 1/ {{piq) {p2<l)\ give infrared singularities (we will 
omit the terms which give finite contributions to the to the decay rate), i.e., 

2 



\M 



2\lR 



aGpVttVZfBC 



'10 



2V2 



TC 



mArca 



+ 16m^ 



p2 



+ 8m 



^ {P2qy 



+ 16m^ 



{piq) , {P2q) 



+ 



(A.2) 



p 



^ipiq)iP2q) ipiq)iP2q) 

Let us consider the first term in the bracket 



d^pi d^p2 d^"- ^^g p'^ 



2Ei 2E2 2 \q\ (27r)"-' {p^qY {2ti 
d^pi 



S{Pb -P1-P2- q) 



p^ 



2E, 2 \q\ (27r)"-^ (p.qY (27r)' 



6 {ml - 2mBEi - 2mBq° + 2Eiq {I - (3z) 



(A.3) 
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where f3 = \p[ \ /Ei, and we choose the first axis along the direction of pi. The integrand 
depends only on the angle 6i = 6, so that we can immediately perform integration over the 
other angles. Integration over all of the angular variables leads to the following result 



1), 



2^n/2-l 



2\ n/2-2 



r(n/2 

where z — cos9. Using (A. 4), and performing integration over radial d \q\ we obtain 



(A.4) 



7 = 



1 



1 27r"/2-i 



2 (27r)"-' (27r)' T (n/2 - 1) Je^-ae AE_ 



47r/? dEi 



dz ■ 



[l-z 



[l-Pz) 



( 2 i^B - '^rnsEi 



, n— 5 



[2ms - 2Ei (1 - 13 z^ 



in— 4 



— 2mB 



{rn?^ — 2mBEi 



(A.5) 



[2mB - 2Ei (1 - 



in— 3 



Second term gives finite contribution and therefore we will omit it. In our case Em = ms/^. 

Introducing next a new variable t — E^ — Ei and expanding all terms in Taylor series 
about El = Ejn, we get 



/2-1)/ 



2n+2^n/2+l p (^/2 - 1) 



dz 



m 
m 



dt- 



(A.6) 



2E^-^{l + (3mz) 



n—A 



After integrating over t we have 

1 



1 



^ 

2n+2^n/2+l p (^/2 - 1) 2 
1 



X 



n - 4 



+ lnAE 



1 (1-2 
dz 



2Nn/2-2 



(A.7) 



n-4 



The last step is to expand (A.7) in a Laurent series about n = 4 and perform trivial 
integrations over dx (retaining only l/(n — 4) and InAii^ terms). Calculating in a similar 
manner all the other terms, we get (we retain only infrared divergent terms and those that 
are proportional to InAE') 



IR 



aCpVt^V^jBCio 



n-4 



2V2 7r 



+ lnAE 



4' 2mB 
8 + 4(l + /?y 



(A.8) 



Now let us calculate O (a) virtual corrections to the B ^ t^t decay. The matrix 
element for the B t~^t~ decay with virtual corrections can be represented as 



M=Mo{l + 47raK + {Z^ - 1)} , 



(A.9) 
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where 



aG p 

Mo = ■^-^fBCioVtbVtl2mr t-^^t , 

is the matrix element for B t'^t" decay without virtual corrections, K denotes O (a) 
corrections due to the photon exchange (vertex) to the Aio, and the last term corresponds 
to the wave function renormalization. Note that all the calculations were performed in the 
Landau gauge. The matrix element square with summation over spins of the final particles 
is given as 

\Mf = \Mof {l + 47rQ;2Re(K) + 2 (Z^ - 1) } . (A.IO) 
In the Landau gauge, fermion wave function renormalization constant is given as 



Zm. — 1 — — 



6 



An Vn - 4 



(A.ll) 



After standard calculation for the infrared singular part for the virtual corrections, we 
finally get 



riR = 



aGFVtkV*fBC,o 



2 

m^mB- 



47r2 



X 



n - 4 



In 



+ 



(A.12) 



+ 



Prom Eqs. (A. 8) and (A.12) we see that the infrared singular terms from Bremstrahlung 
and vertex corrections exactly cancel each other, and the decay rate is infrared-free. 



12 



References 

[1] CLEO Collaboration, R. Ammar et al, Phys. Rev. Lett. 71 (1993) 674; 
CLEO Collaboration, M. S. Alam et al, Phys. Rev.Lett. 74 (1995) 2885. 

[2] A. J. Buras, M. K. Harlander, in Heavy Flavors, edited by A. J. Buras and 
M. Lindner, (World Scientific, Singapore, 1995), pp 58-201. 

[3] Z. Ligeti and M. Wise, Phys. Rev. D53 (1996) 4937. 

[4] G. Buchalla, and A. J. Buras Nucl. Phys. B400 (1993) 225. 

[5] T. M. Aliev, A. Ozpineci, and M.Savci, Phys. Rev. D55 (1997) 7059. 

[6] G. Buchalla, A. J. Buras and M. E. Lautenbacher, Rev. Mod. Phys. 68 (1996) 1125; 
A. J. Buras and M. Miinz, Phys. Rev. D52 (1995) 186. 

[7] M. Misiak, Nucl. Phys. B398 (1993) 23; Erratum: ibid. B439 (1995) 461. 

[8] G. Eilam, I. Halperin and R. R. Mendel Phys. Lett. B361 (1995) 137. 
T. M. Aliev, A. Ozpineci and M. Savci, ibid. B393 (1997) 143. 

[9] Particle Data Group, R. M. Barnett, Phys. Rev. D54 (1996) 1. 

[10] C. S. Lim, T. Morozumi and A. I. Sanda, Phys. Lett. B218 (1989) 343; 

N. G. Deshpande, J. Trampetic and K. Ponose, Phys. Rev. D39 (1989) 1461; 

P. J. O'Donnell and H. K. Tung, Phys. Rev. D43 (1991) 2067; 

P. J. O'Donnell, M. Sutherland and H. K. K. Tung, ibid. D46 (1992) 4091; 

A. I. Vainshtcin, V. 1. Zakharov, L. B. Okun and M. A. Shifman, 

Sov. J. Nucl. Phys. 24 (1976) 427. 



[11] L. Rcina, G. Riccardi, and A. Soni, Phys. Lett. B396 (1997) 231; 
Prep, hep-ph/9706253 (1997). 

[12] W. J. Marciano, and A. Sirlin, Nucl. Phys. B88 (1975) 86. 



13 



